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ABSTRACT 

We use N-body simulations of dark matter haloes in cold dark matter (CDM) and 
a large set of different warm dark matter (WDM) cosmologies to demonstrate that the 
spherically averaged density profile of dark matter haloes has a shape that depends 
on the power spectrum of matter perturbations. Density profiles are steeper in WDM 
but become shallower at r < Virialization isotropizes the velocity dispersion 

in the inner regions of the halo but does not erase the memory of the initial conditions 
in phase space. The location of the observed deviations from CDM in the density 
profile and in phase space can be directly related to the ratio between the halo mass 
and the filtering mass and are most evident in small mass haloes, even for a 34 keV 
thermal relic WDM. The rearrangement of mass within the haloes supports analytic 
models of halo structure that include angular momentum. We also find evidence of 
a dependence of the slope of the inner density profile in CDM cosmologies on the 
halo mass with more massive haloes exhibiting steeper profiles, in agreement with the 
model predictions and with previous simulation results. Our work complements recent 
studies of microhaloes near the filtering scale in CDM and strongly argue against a 
universal shape for the density profile. 
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1 INTRODUCTION 

The seminal work of I Navarro et al.] (1 19961 ) found that the 
density structure of relaxed dark matter haloes are well rep¬ 
resented by what has become known as the NFW profile: 

{r/rs){l + r/rs)^’ 

where p{r) is the density in a spherical shell at distance 
r from the halo centre. By scaling the free parameters Vs 
and ps, which define a characteristic length and density, the 
NFW profile can describe dark matter haloes from dwarf 
galaxy to cluster scales. Furthermore, it was found the NFW 
profile was valid for haloes regardless not only of mass but 
also the power spectrum of initial density fluctuations and 
values of cosmological parameters, establishing that density 
profiles are u niversal in form ind ependent of the cosmolog¬ 
ical context (iNavarro et al.l [19971 ). Another universal prop¬ 
erty was found in the coarse-grained phase-space density 
profile, Q = pja^ , where a is the ve l ocity dispersion of sim¬ 
ulation particles. iTavlor fc Navarr^ (1200 ll ) discovered Q has 
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a remarkably simple form of a power-law, (3 oc r'*', with 
7-1.9. 

It is useful to recast the free parameters of the NFW 
profile in terms of a halo mass and concentration. For relaxed 
haloes a radius can be defined in which the material has 
reached virial equilibrium: 

4:71 Q 

Afvir = A(z)pc('^)R'i;ir5 (2) 

where Rvir is the virial radius enclosing a density A times 
the critical density, pc, at redshift z. The characteristic 
radius Vs can be recast as the concentration parameter, 
Cvir = Rvir It's- Much effort has gone into understanding 
the relationship between Cvir and Mvir as well as the de¬ 
pendencies on the background c osmology and the evolution 
with redshift (IPrada et al.ll2012 ). This is necessary for pre¬ 
dicting the properties of luminous galaxies that reside in 
the dark matter haloes and for using galaxy observations 
as probes of the cold dark matter (CDM) paradigm. The 
concentration was found to correlate with mass such that 
small er mass haloes are more concentrated (INavarro et al.l 
[i 993). This was understood as a consequence of the ear¬ 
lier formation epoch of small mass haloes in the bottom-up 
structure formation of CDM. Since small haloes collapse ear- 
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lier, their inner regions reflect the higher universal density 
of matter at earlier times. Changing the cosmological pa¬ 
rameters or the power spectrum changes the halo formation 
epoch and affects the concentra tions but does not affect th e 
shape of the universal profile (IPolisenskv fc Ricottil [20l3l . 
This interpretation is consistent with simulations of hot and 
warm dark matter which found haloes with masses below 
the truncation scale form later and hav e lower concentra¬ 
tions than CDM halo e s of similar size jiAvila-Reese et al.l 
I 2 OOII: iBode et al.l[200l]: iKnebe et aLll2002lh 

Much effort has also gone into understanding the phys¬ 
ical processes that produce the NFW profile. There are 
two main approaches to analytically modeling the den- 
sity profile: smooth accretion b ased on spherical infall 
dCunn fc Gotti Il972l : ICottl Il975ll . and hierarchical merg- 


ing following Press-Schechter formalism (IPress & Schechtei 

1974 

2003 

: IPeebled Il974l: iLacev & Cold Il993l: Manriaue et al. 

). Both approaches have been successful at producing 


the universal profiles. This has been explained as a result of 
the process of virialization. If virialization erases all infor¬ 
mation about the past merging history of the halo then it 
does not matter if the mass accretion is modeled as clumpy 
or smooth. However, a consensus has not emerged on the 
dominant processes occurring during virialization or if the 
virialization process erases all memory of the initial condi¬ 
tions. 


The early stages of halo formation are marked by rapid 
accretion and mergers making it natural to consider vio¬ 
lent relaxation as the dominant mechanism determining the 
dar k matter pro files in the fluctuating gravitational poten¬ 
tial (IWhitell996l 'l. Violent relaxation was origin ally proposed 
to ex plain the structure of elliptical galaxies dLvnden-Belll 
[1963) where estimates of star-star encounters would not es¬ 
tablish equilibrium in a Hubble time. The relaxation time of 
a forming halo is re lated to the r a te of chang e of the grav¬ 
itation al potential. lAustin et al.l |2003) and I Barnes et al.l 
d2006tl argue the universal nature of Q{r) results from vi¬ 
olent relaxation. 

The works of lWechsler et af] d2002l '): [Zhao et ^ d2003l . 
l2009l i have shown there are two main eras of halo growth, a 
fast accretion phase and a slow phase. The fast growth phase 
in CDM is dominated by mergers of objects with similar 
mass in contrast to the slow growth phase characterized by 
quiescent accretion and minor mergers. The inner halo is 
set at the end of the fast era with the slow growth phase 
having little impact on the inner structure and gravitational 
potential well, leading to an inside-out growth of haloes. 
These studies And violent relaxation is only important in 
forming the inner profile with the outer profile determined 
by secondary infall during the slow growth phase. 


The NFW profile is characterized by a logarithmic 
slope, a = dlog p/dlogr, that rolls from an asymptotic 
value a = —3 at large radii to a = —1 in the inner 
halo. The value of the inner slope has been a matter of 
controversy. The first concerne d the value of the asymp¬ 
totic slope dMoore et al.l Il999li . As the number of parti¬ 
cles in simulations have increased it has become evident 
the density profiles do not approach an asymptotic value 
but continue to roll slowly with radius ([Navarro et al.ll2004l : 
iDiemand et al.ir2004l : lGr aharn_eir^ [gOOgh and are better de¬ 
scribed by Einasto profiles ( Einastol 1965h . However, this has 


not changed the conclusion that all information about the 
formation history is lost in the virialization process. 

The second controversy is a dependence of the inner 
profile on the halo mass. Many models have been con¬ 
structed that explain the emergence of the universal profile 
as a consequence of repeated mergers (ISv er fc White! Il9^: 
Nusser fc Shethlll999l : [^bramanian et S' 2000l : Dekel et ah! 


2 OO 3 II . Although they differ in the details, the relevant phys¬ 


ical processes determining the halo properties are the tidal 
stripping of material from accreting subhaloes, dynamical 
friction, and tidal compression transferring energy from the 
satellites to the halo particles and the decaying of satellite 
orbits to the halo centre. These models predict a dependence 
of the inner density profile slope on the slope of the power 
spectrum at the scale of the halo, P{k) oc fc". The steeper 
spectrum characteristic of dwarf-scales is predicted to pro¬ 
duce softer cores than for galactic and cluster-scale h aloes. 
Independent of the merger models, iDel Popold ll2010l l ques¬ 
tions the universality of both the density and the Q profile 
and concludes both should depend on mass. His spherical 
infall models with angular momentum show a steepening of 
the inner density profile with increasing halo mass, although 
to a lesser extent than the merger models. The central ques¬ 
tion of these studies is, do haloes in equilibrium retain any 
memory of the initial conditions and mass function of ac¬ 
creting satellites they are built from or is all information 
lost in the virialization process? 

iRicottil (l2003f) ran CDM simulations of the same real¬ 
ization of the density field in boxes of varying side length to 
compare the profiles at different mass scales. He examined 
the average profiles when the box structures showed similar 
clustering and the most massive haloes were composed of 
the same number of particles. He found a systematic depen¬ 
dence of the inner slope on halo mass with dwarf-scale haloes 
having softer cores than galactic and cluster-scale haloes 


in agreement with the predictions of Subrarnanianet^ 


^re 

3 ). 


(200d) . These results were reinforced iniRicotti fc Wilkinson! 


2004) ; iRicotti et al.l (l2007ll . Ijing fc Sutd ll2000ll also saw a 


dependence of inner slope in their simulations of haloes at 
galactic and cluster scales. 

Another way of testing the importance of substructure 
is by introducing a truncation in the power spectrum as 
in hot and warm dark matter (WDM) cosmologies where 
substructure is suppressed below the particle free-streaming 
scale and haloes form by monolithic collapse. Many in- 
vesti gations using these c o smologies have be e n conducted 


vesti gations using tnese c o smologies nave be e n conducted 
(e.g. ]WarigL_fcJwiiE3 2009I. Huss_et^ _ 1999bl. Mooreet^D 


1999l. lColm et al.ll2000l. iBusha et al.ll2007l . iBode et~ I l200ll . 


Polisensky fc Ricottil 201ll '). These works find haloes that 


form below the truncation mass have lower concentrations 
consistent with their later formation epochs but the profiles 
ar e well described by the NFW form. This is in contrast 
to IColfn et al.l (I2008II who find the inner profiles are sys¬ 
tematically steeper in their WDM simulations of Galaxy¬ 
sized haloes. Recent studies of CDM microhaloes confirm 
halo profiles near the streaming scale are steeper when 
a truncation is introduced into the power spectrum than 
witho ut, however the origin of the steepening remains un¬ 
clear (lls hivama et al.l 120101 : lAnderhalden fc Diemandll2013l : 
Ilshivamall20l4l . Independent models of the density profiles 
of WDM haloes predict not a steepening but a flattening 
of the inner profile due to the truncated power spectrum 
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dWilliams et al.|[2004l : IVifias et aDl2012l ~). It is important to 
stress the flattening in these models is due to the truncated 
power spectrum alone, not to the random thermal motions 
of WDM particles which were not included in the models. 

In this work we employ N-body simulations of halo for¬ 
mation in CDM and WDM cosmologies to explore the effects 
of the power spect rum on halo st ructure and dynamics. We 
use the method of I Rico ttl (l2003h of scaling the simulation 
volume to change the mass scale (§ El). However, we do not 
study a statistical sample of haloes but focus on a handful 
of halos with about a factor of 100 greater mass resolution 
than in iR.icottil (l200.ll 'l. The goal is not to rigorously test 
any particular halo model but simply to look for evidence 
the halo retains memory of the initial power spectrum. This 
evidence is expected to manifest itself as trends in the halo 
profiles as the mass and truncation scales change (§ ISl). We 
examine if our results are typical of a larger halo population 
(§11 and investigate the physical origin of the results (§ El). 
A discussion and comparison to previous published works is 
presented in § [6] and the summary in § [T] 

2 NUMERICAL SIMULATIONS 
2.1 Cosmological Models 

WDM particles are relativistic in the early universe and free- 
stream out of overdense regions before the adiabatic expan¬ 
sion of the universe reduces the particles to subrelativistic 
velocities. WDM thus damps density perturbations below a 
characteristic scale that depends on the particle mass and 
acts as a filter on the power spectrum of density perturba¬ 
tions. The power spectra for WDM cosmologies is related to 
that for CDM by 

Pw{k) = PcT^, (3) 

where Tw is the WDM transfer function. The transfer func¬ 
tion given bv lBode et al.l ((2003) is used for dark matter par¬ 
ticles that were coupled to the relativistic cosmic plasma at 
early times and achieved thermal equilibr ium prior to the 
time o f their decoupling. The formula of lEisenstein fc Hul 
lll998tl is adopted for the CDM power spectrum. 

We define the WDM filtering mass as in 
ISommer-Larsen fc Dolg^ ll200ll '). 

_q 

^ t^mpck t , (4) 

3 

where pc is the critical density and fc/ is a characteristic free- 
streaming wave number defined where Tw = 0.5. For consis¬ 
tency with ISommer-Larsen fc Dolgov! ll200lh we also define 
the free-streaming, or Altering length as Rf = OAQkJ^. 

We adopted values for cosmological p arameters from 
the Bolshoi simulation (|KlvDin et al.llioill 'l. (Dm, Ha, H(,, 
h, as, ria) = (0.27, 0.73, 0.0469, 0.7, 0.82, 0.95), which were 
chosen to be within la of WMAP5, WMAP7, and consistent 
with the results of supernovae, and X-ray cluster surveys. 
These parameters are also within 1.7 ct of WMAP9 and 2.2 ct 
of Planckl. We use a variety of WDM models for thermal 
relics in the range 0.75 — 66 keV. 

Since our focus is to examine the effects of the power 
spectrum on halo structure, the initial conditions include 
particle velocities due to the gravitational potential using 
the Zeldovich approximation but random thermal velocities 


appropriate for WDM have not been added to the simula¬ 
tion particles. For the WDM cosmologies adopted here the 
effects of thermal velocities are expected to be small; this is 
discussed further in Section El 


2.2 Software 

The simulations were conducte d with the A-b ody cosmolog¬ 
ical simulation code GADGET -2 (ISDringelll2005l l with gravita¬ 
tional physics only and init ial conditions gener ated with the 
GRAFIC 2 software package llBertschingeiluOOlll . We produce 
a single realization of the density field but vary the power 
spectrum of fluctuations appropriate for CDM and WDM 
cosmologies. 

The AMIGA’S Halo Finder (ahf) software 
llKnolImann fc Knebd l2009ll was used to identify and 
characterize all gravitationally bound haloes composed of 
at least 50 particles after iteratively removing unbound 
particles. The virial mass of a halo is defined in Equation El 
Since the simulations are confined to high redshifts {z > 3) 
the universe is matter dominated at all epochs and we 
adopt the virial condition for an Einstein-de Sitter cos¬ 
mology, A{z) — 178. The MergerTree tool in AHF was used 
to construct merger trees, identify halo progenitors at all 
times, and for identifying haloes across cosmologies. 

AHF calcul ates the convergen ce radius according to 
the criterion of IPower et al.l (l2003l l and is generally about 
10 softening lengths, enclosing ~ 2000 particles at r ~ 
0.00677^ir. We tested this by running low resolution sim¬ 
ulations and found that the profiles are actually converged 
to about 5 — 6 softening lengths, enclosing ~ 200 particles at 
r ~ 0.003i?„ir. The convergence radius given by AHF may 
be overly conservative for the simulations but this has no 
impact on the results. When examining the halo profiles we 
adopt the convention of plotting r > 6e but indicating in 
bold where the profiles satisfy the criterion of IPower et al.l 


2.3 Simulations 

We simulated a small cubic box with a comoving side length 
of 3.3 Mpc from z = 79 to z = 8 with 512® particles and 
mass resolution ~ 10 ^Mq. A halo of mass ~ 2x 10 ®Mq that 
appeared to have an early formation epoch and relaxed to 
virial equilibrium at scale factor a = 0.1 {z = — 1) was 

chosen for resimulation using a zoom technique. We refer 
to this as “Halo A.” A volume of higher mass resolution 
was generated in the initial conditions covering the initial 
volume of particles that end within three virial radii of Halo 
A. We ran high resolution simulations with a mass resolution 
of 18.7 Mq in CDM and multiple WDM cosmologies in the 
range 4—66 keV. We adopted a force softening length of 8 pc, 
held constant in comoving units. To test the convergence of 
our results we also ran low resolution simulations with the 
mass resolution reduced by a factor of 8 in CDM and 6 keV 
WDM. We further tested the dependence of our results on 
the initial conditions by running low resolution tests in CDM 
and 6 keV starting from 2 = 120. 

We ran additional simulations with the box size in¬ 
creased to medium and large side lengths of 7.0 and 22.4 Mpc 
to increase the mass scale a factor of 10 and 320, respectively. 
The force softening lengths were also scaled with the box size 
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redshift 



0.05 0.1 0.15 


scale factor 

Figure 1. Mass growth of Halo A in the small, medium, and 
large box simulations (top to bottom) in CDM and select WDM 
cosmologies. The circles show the CDM halo at the normalization 
times when the halo has grown to ~ 10^ particles at the three 
mass scales. 

to 17 pc and 55 pc, respectively. CDM and WDM cosmolo¬ 
gies ranging from 2 — 5 keV were run for the medium mass 
scale while CDM and 0.75 — 2 keV WDM were run for the 
large mass scale. To compare Halo A across mass scales we 
define “normalization times” as the epochs when the CDM 
haloes have grown to encompass the same number of par¬ 
ticles within the virial radius, N ~ 10^, as the small mass 
scale at a = 0.1. This occurred at a = 0.116 and a = 0.155 
for the medium and large scales, respectively. Figure[T] shows 
the growth of Halo A for the three mass scales. Halo forma¬ 
tion Is delayed in the WDM cosmologies. However, once it 
begins it grows quickly until it catches up with the CDM 
halo, after which it evolves at a similar rate. The circles 
show the normalization times when the halo has entered the 
slow-growth phase and has reached the same size in CDM 
at all scales. At the normalization times the halo masses are 
approximately 2 x 10* Mq, 2 x 10® Mq, and 6 x 10^® Mq for 
the small, medium, and large mass scales, respectively. 

Table [T] gives a summary of the simulations conducted 
in this work. Listed in the table are the WDM Hltering mass, 
the Hltering length, and the number of simulation particles 
sampling the filtering mass, Nf. For the WDM cosmologies 
the box side length Li,ox is given in units of the filtering 
length. This is a convenient way to show in which simu¬ 
lations the effects of the truncated power spectrum will be 
similar across mass scales. The 2 keV simulation in the large 
box is expected to be similar to the 5 keV medium box and 
10 keV small box. Likewise, the 1.1 keV large box will be 


Table 1. Properties of simulations. 


Cosmo 

Mf 

Rf 

Nf 

Lb. 

OX 


[Mq] 

[kpc] 

[xlO'^] 



(1) 

(2) 

(3) 

(4) 

(5) 


Small Box: mres = 

: 18.7 Mq 



CDM 

- 

_ 

_ 

3270 

kpc 

66 keV 

1.90 X 10^ 

0.7 

lO-'^ 

4485 

Rf 

34 keV 

1.90 X 10® 

1.6 

0.001 

2082 

Rf 

21 keV 

9.50 X 10® 

2.7 

0.005 

1217 

Rf 

15 keV 

3.19 X 10® 

4.0 

0.02 

812 

Rf 

10 keV 

1.29 X 10’^ 

6.4 

0.07 

510 

Rf 

7 keV 

4.43 X 10'^ 

9.7 

0.2 

338 

Rf 

6 keV 

7.54 X 10'^ 

11.5 

0.4 

283 

Rf 

5 keV 

1.41 X 10® 

14.2 

0.8 

230 

Rf 

4.5 keV 

2.03 X 10® 

16.1 

1.1 

203 

Rf 

4 keV 

3.05 X 10® 

18.4 

1.6 

178 

Rf 


Medium Box: 

TTifes 

= 187 Mq 



CDM 

- 

_ 

_ 

7049 

kpc 

5 keV 

1.41 X 10® 

14.2 

0.08 

495 

Rf 

3 keV 

8.24 X 10® 

25.6 

0.4 

275 

Rf 

2 keV 

3.34 X 10® 

40.8 

1.8 

173 

Rf 


Large Box: m 

res — 

6, 010 Mq 



CDM 

_ 

_ 

_ 

22394 kpc 

2 keV 

3.34 X 10® 

40.8 

0.06 

548 

Rf 

1.1 keV 

2.63 X lOi® 

81.2 

0.4 

276 

Rf 

0.75 keV 

9.84 X lO^® 

126.2 

1.6 

177 

Rf 


similar to the 3 keV medium box and 6 keV small box. The 
0.75 keV large box will be similar to the 2 keV medium box 
and 4 keV small box. These simulations are colour-coded 
in Figure [T] according to their ratios of box scale to WDM 
filtering scale. Another way of characterizing the similarity 
of these simulations is by the ratio of filtering mass to halo 
virial mass. This ratio is equivalent to Nf expressed in units 
of 10^, as listed in column (4) of Table[T] For the similar cos¬ 
mologies given above, the filtering masses are approximately 
7%, 40%, and 170% of the halo masses at the normalization 
times. 

Table [2] summarizes the properties of Halo A at the nor¬ 
malization times. An examination of Figure [T] shows Halo A 
has not suffered a recent major merger and is in the slow 
growth phase in all cosmologies at the normalization times. 
However, a more rigourous examination of the halo relax¬ 
ation state is desirable. Differences from a universal profile 
are seen in unre laxed haloe s and haloes with large amounts 
of substructure (Ijln aH). Additionally, the Inner slope of 
the density profile is sensitive to the location of the halo 
centre. An artificial flattening of the profile could be pro¬ 
duced by an ambiguously defined centre due to a recently 
arrived subhalo at the core, for example. We performed a 
qualitative visual examination that the halo centres deter¬ 
mined by AHF correspond to the density peak of particles 
and we examined quantitative measures of the relaxation. 
Studies with large samples of haloes have identi fied several 
metr i cs for separating haloes b y relaxation state dNeto et al.l 
I 2 OO 7 I : [Maccio et ^l2007l . l2008l ): Xoff, the offset between the 
halo centre and centre of mass of particles within Rvir', the 
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Table 2. Properties of Halo A at the normalization times in the 
high resolution simulations. 


Cosmo 

ALyir 

A' 


2K 1 

|/ 7 | ^ 

fsub 


[10®Mo] 

[10-2] 



(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

Halo A - Small Box 

CDM 

1.868 

4.21 

0.06 

0.41 

0.08 

66 keV 

1.847 

3.84 

0.05 

0.40 

0.07 

34 keV 

1.863 

3.89 

0.05 

0.40 

0.06 

21 keV 

1.882 

3.99 

0.06 

0.39 

0.05 

15 keV 

1.888 

4.24 

0.05 

0.42 

0.05 

10 keV 

1.887 

4.63 

0.08 

0.42 

0.05 

7 keV 

1.811 

4.88 

0.08 

0.40 

0.04 

6 keV 

1.713 

4.80 

0.07 

0.40 

0.03 

5 keV 

1.507 

4.30 

0.06 

0.39 

0.03 

4.5 keV 

1.330 

3.93 

0.06 

0.39 

0.03 

4 keV 

1.074 

3.16 

0.06 

0.41 

0.05 

Halo A - Medium Box 

CDM 

18.533 

4.09 

0.06 

0.39 

0.08 

5 keV 

18.894 

4.63 

0.06 

0.41 

0.05 

3 keV 

16.912 

4.76 

0.07 

0.40 

0.04 

2 keV 

10.126 

3.10 

0.05 

0.40 

0.05 

Halo A - Large Box 

CDM 

600.497 

3.34 

0.04 

0.39 

0.13 

2 keV 

604.681 

4.07 

0.08 

0.40 

0.05 

1 keV 

547.367 

4.44 

0.07 

0.39 

0.03 

0.75 keV 

351.273 

3.48 

0.04 

0.38 

0.05 


virial ratio 2K/\U\ — 1; the mass fraction bound in subhaloes 
fsuh\ and the spin pa rameter, A^ that ch aracterizes the halo 
angular momentum (jBullock et alJlioOll ): 


y/2MvirV vir K vir 

where J is the total angular momentum of all particles 
within Rvir and Vvir is the circular velocity at Rvir, = 
GM/R. These metrics are listed in Table [2] The general 
conditions for a relaxed halo are; A' < 0.1, Xpff < G.QlRyjr , 
2KI \U\-1 < 0.35, and fsub < 0.1 (|Ludlow et al.ll20lll2013l. 
I 2 OI 4 I 1 . Halo A largely satisfies these criteria with the great¬ 
est discrepancy being a slightly larger virial ratio of ~ 0.4, 
however in Section [4] we examine the stability of the profiles 
and find our results are not due to the relaxation state of 
the halo or transient accretion events. 

We note a curious increase in fs^b for the warmest cos¬ 
mologies in all three box sizes. We examined this further by 
applying varying subhalo mass cuts and found faub is dom¬ 
inated by subhaloes with less than 5000 particles in these 
cosmologies. Artifical haloes of this size are known to form 
along filaments in tru ncated power spectrum cosmologies 
dWang fc Whitell2007l ~) and are likely contaminating faub- 


3 RESULTS I - NON-UNIVERSALITY OF 
PROFILES 

We begin by examining the effects of the WDM power spec¬ 
tra on the density structure of Halo A in the three boxes and 
thus the three mass scales of the halo. We then examine the 
kinematics and conclude by checking the convergence. In 
Section |4] we examine a larger halo sample to check if the 
results of Halo A are typical for haloes in general and we 
show the features in the profiles are dynamically stable. 


3.1 Density Structure 

The left panel of Figure [5] shows the spherically averaged 
density profiles of Halo A at a = 0.1 in all cosmologies for 
the small mass simulations, M — 2x 10® Mq. The profiles are 
plotted wit h solid lines where they satisfy the convergence 
criterion of I Power et all ll2003ll and the inner profiles are ex¬ 
tended to six force softening lengths with dotted lines. Excel¬ 
lent agreement is seen between the high and low resolution 
CDM profiles. No significant differences from CDM are seen 
for a filtering mass 10“"^ of the halo mass. As the filtering 
mass increases the inner profile at r < O.lRuir steepens and 
the density increases. The location where the WDM density 
begins to increase shows a correlation with filtering scale, 
moving to larger radii as the cosmology gets warmer and the 
filtering scale larger. In the range 0.1 — 0.4Ruir the WDM 
densities fall below that of CDM while the densities agree 
with CDM for r > OARvir except in the warmest cosmolo¬ 
gies where the halo outskirts are still growing by secondary 
infall. These features are more pronounced in the cumula¬ 
tive mass profiles shown in the right panel of Figure [2l The 
enclosed mass is equivalent in CDM and WDM at > 0.5Rvir 
indicating it is the mass in shells at 0.1 — O.^Ruir that has 
been displaced to smaller radii in the WDM simulations. 

To compare the mass profiles of Halo A across the three 
mass scales we plot in the left panel of Figure [3] the profiles 
of enclosed number of simulation particles and normalize the 
radial coordinates by the CDM virial radius in each box. In¬ 
terestingly, variations are seen at r < O.lRuir in the CDM 
haloes as a function of the halo mass in contrast to the WDM 
simulations where the profiles are nearly identical across the 
explored halo mass range. The enclosed mass in the CDM 
inner halo becomes greater as the halo mass increases but 
when small scale structures are erased, as in the WDM sim¬ 
ulations, the mass profiles are insensitive to the halo mass. 
Angular mo mentum sets the s hape of the inner profile in 
the models of lPel Popolol l|2nn9h where more massive haloes 
are predicted to have less angular momentum resulting in 
steeper profiles. It can be seen from Table [2]that the spin of 
the CDM halo decreases as the mass scale increases, consis¬ 
tent with this idea. 

The CDM halo spin parameter is 26% higher at the 
small mass scale compared to the large while the WDM 
haloes vary by ^ 10% which may be why the WDM pro¬ 
files are very similar. However, the WDM spin parameters 
are generally higher than the CDM halo at all scales yet 
they have steeper profiles than CDM so this is not the en¬ 
tire answer. One difficulty is the spin parameter includes all 
particles within the virial radius while the greatest differ¬ 
ences between profiles are seen in the inner halo. This is 
examined further in Section (5] 
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Figure 2. (Left) Density profiles of the small mass simulations of Halo A at the normalization time, a = 0.1. The dashed line gives 
the asymptotic slope of the NFW profile at small radii. (Right) Cumulative mass profiles of the small mass simulations of Halo A. The 
density profiles have been multiplied by P and the mass profiles by r~^ to reduce the dynamic range. The radial coordinates have 
been normalized to the virial radius in CDM and are plotted to six softening lengths (fie). The WDM profiles have been grouped and 
are plotted against the CDM profile for clarity. All profiles are plotted with solid lines where they satisfy the convergence criterion of 
I Power et al.l ll200.lh . As the filtering mass increases clear deviations from a universal shape are seen in the WDM profiles caused by the 
displacement of mass from intermediate radii to the core in the WDM simulations. 


It is important to emphasize that the differences be¬ 
tween WDM and CDM profiles diminish as the halo mass 
increases due to the steepening of the CDM profile. This ob¬ 
servation may explain why previous works have not clearly 
identified the prominent features and trends in the profile 
shapes found in this work and the rec ent simulations of mi¬ 
crohaloes near the CDM filtering scale llls hivania et al.l[201^ : 
lAnderhalden fc Diemandll2013l : Tlshivamall2014l ). 

In the right panel of Figure [3] the logarithmic slope of 
the density profiles are compared across mass scales for a 
given ratio Mf/M. The large mass CDM halo profile is 
steeper than the medium and small haloes for r < O.SRvir 
and reaches the NFW value of —1 at a smaller radius (given 
by the vertical short grey lines). Ftowever, the differences are 
le ss than predicted by the m odel of undigested subhalo cores 
of ISubramanian et al.l (120001 ) but in agreement with the pre¬ 
dictio ns of the angular momentum models of iDel Popold 
II 2 OIOI) . 

Unlike the CDM profiles the slopes in the WDM cos¬ 
mologies are nearly identical across mass scales. For Mf < 
M the slopes tend toward plateaus of constant value when 
moving to smaller radii. These plateaus also depend on 
the filtering scale, ranging from r = 0.02 — O.lRvir for 
Mf ~ 0.07M but grow steeper and move outward to r = 
0.04 — O.SRvir for Mf ^ QAM . These ten dencies w e re als o 
seen in the CDM microhalo simulations of llshivamal ll2014h . 
They fit a power law function to the relation between halo 
mass and inner slope a for their sample: 

a = 0.123 log(M„i^/M/) - 1.461, ( 6 ) 

with a scatter of 20%. We plot the values given by this equa¬ 
tion as horizontal arrows in the right panel of Figure 0 and 
find agreement with the slope plateaus to 10 %. 

Our simulations also show inte r esting new features in 
regimes not explored in llshivam^ (|2014|) . Continuing to 
smaller radii the slopes do not remain at the asymptotic 
values of the plateaus but soften. The slopes remain steeper 


than CDM for r < O.lRvir and achieve —1 at smaller radii, 
although this scale moves outward as the filtering scale 
gets larger. The inner profiles quickly become softer than 
CDM at r ^ Q.OlRvir- The inner profile s are shallower 
in WD M i n agreement with the models of IWilliams et aP 
ll2004l) and IVifias et al.l ll2012l) . For Mf > M the slope pro¬ 
files do not form a plateau but remain steeper than CDM 
for r > 0.01i?„ir. 


3.2 Internal Kinematics 

Figure |4] shows the profiles of P for the small mass simula¬ 
tions of Halo A, where a is the local 3D velocity dispersion. 
Similar to the density profiles, the dispersions are greater 
in the inner WDM haloes compared to CDM and show a 
correlation with the filtering scale, growing larger and ex¬ 
tending to greater radii as the cosmology becomes warmer 
and the filtering scale increases. This can be understood as 
a consequence of the increased mass in the WDM cores. As 
the mass in the core grows the dispersion must get larger to 
stay in virial equilibrium against the deeper potential well. 

To examine the phase space density profiles of Halo A 
we adopt 7 = —1.875 and fit the profiles of each CDM 
halo to the form, Qfu = AM. For illustrative purposes 
we show Q normalized to the power law fit, Q = QjQfu, 
to emphasize deviations from a simple power law. We cal¬ 
culate Q for the WDM simulations using the CDM Qfu. 
The left panel of Figure [S] shows the deviations from power 
law for Halo A in common cosmologies for all three mass 
scales. A prominent feature is seen in the inner regions of 
the WDM haloes, reaching a maximum deviation before de¬ 
clining. Along the top axis of each plot are ticks marking the 
location of 0.03777/ in each WDM cosmology. This scaling 
was empirically determined but marks the location of the 
peak remarkably well indicating the deviations scale with 
the filtering scale. 








































Fingerprints of the ICs 7 



0.01 0.1 1 0.01 0.1 1 


^/^vir.CDM 



’’/^vir.CDM 


Figure 3. (Left) Comparison of the cumulative mass profiles of Halo A at the normalization times of the three mass scales when the 
haloes have grown to 10^ particles. The small, medium, and large mass haloes are plotted with the solid, dotted, and dashed lines, 
respectively. The profiles are given by the number of enclosed simulation particles and the radial coordinates have been normalized by 
the CDM virial radii. The WDM profiles are plotted against the CDM profiles for comparison. (Right) Comparison of the slope of the 
density profiles of Halo A across the three mass scales at the normalization times. Short grey lines and vertical arrows indicate where the 
log slope is —1 in all cosmologies. The CDM halo is steeper at the large mass scale than the smal l scale. The WD M haloes are generally 
steeper than the CDM but soften at r < 0.01Rt,ir- Horizontal arrows show the fitting formula of llshivamal |[2o3) match the WDM 
plateaus to 10%. 



’'/^vir.CDM 


Figure 4. Velocity dispersion profiles of the small mass simula¬ 
tions of Halo A. WDM simulations have been grouped and plotted 
against the CDM profile for clarity. 


Interestingly, a drop in the inner profile of the CDM 
haloes is seen that becomes more pronounced as the 
halo mass decrea ses. This also agrees with the models of 
iDel PoDold (|20inh that demonstrate a dependence of the Q 
profile on halo mass. 

A useful metric of the particle orbits is the velocity 
anisotropy parameter given by: 


I3{r) = 1 - 


I „2 

2cr2 


(7) 


where ag and are the angular velocity dispersions and 
is the radial velocity dispersion. For purely radial orbits, 
P = 1, while isotropic particle motions give /3 = 0. At the 
halo outskirts /? —>■ 1 where freshly accreted material is still 
falling inward while in the core of a relaxed halo /3 —>■ 0. 


In practice the anisotropy parameter is seldom exactly zero 
in the inner regions since simulated haloes are generally not 
spherically symmetric. 

The right panel of Figure [5] shows the velocity 
anisotropy profiles of Halo A at common filtering scales in 
all simulations. There is a radial bias in the CDM particle 
orbits at r > O.lR^ir while particles inside this scale are 
well isotropized. The WDM anisotropy profiles are gener¬ 
ally similar to CDM, although in the warmest cosmologies 
(bottom row) the haloes have radial bias extending deeper 
into the inner halo than in the other cosmologies. An exam¬ 
ination of the halo shapes revealed a tendency for the haloes 
to become more spherical in the inner and outer regions at 
r < 0.02Rvir and r > Q.2R^ir and less spherical in the inter¬ 
mediate regions as the filtering mass increases. The features 
of the velocity anisotropy may thus not be due to an incom¬ 
plete isotropization of the particle velocities but simply to 
the increased triaxiality seen at these radii. What is clear is 
the lack of any feature at the location of the peak deviations 
in the Q profile making it apparent the physical processes 
that created the increased mass in the WDM cores do not 
leave an imprint on the isotropy of particle velocities after 
virialization. 


3.3 Convergence Tests 

Figure^] shows the profiles of Halo A in the high and low res¬ 
olution simulations and the test simulations initiated from 
a higher redshift. Excellent agreement is seen across reso¬ 
lutions in both CDM and WDM and from the simulations 
started from highe r redshift. The convergence criterion of 
I Power et al.l (l2003l l tested on CDM simulations appears to 
be not only valid but perhaps overly conservative in mea¬ 
suring the convergence radii of our WDM simulations. 

Simulations with truncated power spectra are known to 
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Figure 5. (Left) Deviations from power-law behaviour in the phase space density profiles of the simulations of Halo A. (Right) Velocity 
anisotropy profiles of the simulations of Halo A. The simulations are grouped by mass scale: small, medium, and large from left to right, 
and by relation of filtering mass to halo mass with the cosmology growing warmer from top to bottom. Deviations from power-law are 
seen in the inner WDM haloes that reach a peak at 4% of the filtering length (coloured ticks along the top of each plot). However, no 
corresponding features are seen in the velocity anisotropy. 



0.1 1 10 


Figure 6. Comparison of the high and low resolution small mass 
simulations of Halo A. The high and low resolution CDM simula¬ 
tions are plotted in black and grey, respectively, the high and low 
6 keV simulations in light and dark blue. The dashed lines are the 
low resolution 6 keV and CDM simulations started from a higher 
redshift. Profiles are plotte d with solid lines w here they satisfy 
the convergence criterion of I Power et al.l ll2003h and are extended 
to 3e with dotted lines. Consistent results are seen across sim¬ 
ulations demonstrating the results are not affected by the mass 
resolution or starting redshift. 


produce numerically artificial small mass hal oes along the fil¬ 
amen ts of collapsed density perturbations (|Wang fc Whitel 
|2007^ whose size and separation are dependent on the mass 
resolution. Figure |6] also shows the results are not due to 
these spurious haloes since a dependence on the mass reso¬ 
lution would be expected to reflect on the shape and location 
of the features in the WDM profiles. 


4 RESULTS II - TESTING COSMIC 

VARIANCE 

To explore if the results for Halo A were typical of haloes 
in general we simulated two additional haloes, Halo B and 
Halo C. We first ran a low resolution simulation with a cu¬ 
bic refinement volume of side length 1/4 the box length, 
composed of 512® particles, and centred on Halo A. The 
15 largest haloes in the rehnement volume were examined 
in detail. Two haloes were chosen based on their quiescent 
accretion histories and relaxation metrics for individual res¬ 
imulation at high resolution in an analogous way to Halo 
A but only at large mass scale and only in CDM, 1.1 keV 
and 0.75 keV cosmologies. Haloes B and C grow to about 
the same size as Halo A by the end of the simulations but 
have later formation epochs suggesting Halo A forms from 
a volume with greater initial overdensity. For example, the 
0.75 keV large scale Halo A grows to A = 10® at a = 0.125 
while Haloes B and C don’t reach this size until a = 0.148 
and 0.165, respectively. 

We ran the large scale simulations of Haloes B and C 
to a = 0.25 to test the stability of the profiles over a greater 
time span than the simulations on Halo A presented in Sec¬ 
tion [3] To test the dynamical stability of Halo A we ran an 
extended simulation set at small mass scale for CDM, 6 keV 
and 4 keV WDM. The refinement volume was increased a 
factor of two to allow Halo A to evolve to a = 0.125 while 
preventing contamination by accretion of particles from out¬ 
side the rehnement volume. Snapshots of the particle data 
were saved every 0.001 change in scale factor. We list the 
properties of all three haloes at the end of their simulations 
in Table |3l 

Figure [7] shows the density prohles and slopes for Haloes 
B and C and give excellent agreement with the results from 
Halo A. The slopes of the WDM prohles show the same ten¬ 
dency to steepen, reach a plateau, then soften in the inner 
halo, and show the same dependence on hltering scale. The 
density prohles similarly show the same characteristic fea- 
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Table 3. Properties of the small scale Halo A and large scale 
Haloes B and C at the end of each extended simulation set. 


Cosmo 

A'lviv 

A' 

^off 

2K 1 

W\~^ 

fsub 


[1O®M0] 

[10-2] 

[Rvir] 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


Halo A - Small Box 

- a = 0.125 


A CDM 

3.789 

4.87 

0.07 

0.41 

0.04 

A 6 keV 

4.029 

4.86 

0.04 

0.40 

0.01 

A 4 keV 

3.362 

3.25 

0.04 

0.37 

0.01 

Haloes B & C 

- Large Box - a 

= 0.25 


B CDM 

802.734 

1.96 

0.03 

0.24 

0.09 

B 1.1 keV 

708.649 

1.77 

0.02 

0.21 

0.01 

B 0.75 keV 

565.010 

1.23 

0.04 

0.21 

0.01 

C CDM 

572.440 

2.72 

0.05 

0.21 

0.12 

C 1.1 keV 

451.046 

2.78 

0.06 

0.19 

0.02 

C 0.75 keV 

308.477 

1.45 

0.02 

0.17 

0.01 



^/^vir.CDM 


Figure 7. (Top) Density profiles of Haloes B and C in the large 
box-size simulations. The density profiles have been multiplied by 
to reduce the dynamic range. (Bottom) Slope of the density 
profiles for the same two halos. 

ture of mass displaced from intermediate regions to the core 
in the WDM cosmologies. The main difference is a stronger 
deviation from CDM with a larger overall reduction in den¬ 
sities. The density profiles are everywhere below the cor¬ 
responding CDM halo for the 0.75 keV cosmologies. Halo 
C shows a larger deviation from CDM than Halo B, likely 
related to the later formation epochs of these haloes. 

We show the normalized phase space density profiles, 
Q, in Figure [ 8 ] for all three haloes in their CDM and WDM 
simulations at times when the haloes are composed of > 10 ® 
particles. The values of 7 for Qfu were determined sepa¬ 
rately for each halo by fitting the CDM profiles at late times. 
The CDM profiles of all three haloes are consistent with mi¬ 
nor fluctuations about power laws while the large deviation 
seen in the WDM core of Halo A and its dependence on the 
filtering scale is also seen in the other two haloes. 

The prominent feature of the 6 keV WDM core first 
appears in the profile of small scale Halo A at a = 0.075, 


and at a = 0.115 in the large scale 1.1 keV Haloes B and 
C. After its formation the core feature is stable through the 
end of the simulation, spanning a time range of 420 Myr for 
Halo A and 1520 Myr for B and C. The dynamical time, 
tdyn = \JB?/GM, for the core of Halo A is 6.5 Myr and 
14-18 Myr for the cores of Haloes B and C, demonstrating 
the profiles are stable over at least 65 — 85 dynamical times. 
In the warmer cosmologies the cores form later but are also 
stable through the end of the simulations. 

The mass growth rates, M = d log M/d log a, when the 
core features appear in the WDM profiles are M = 6 — 
10. At earlier times the growth rates are higher and the 
profiles fluctuate wildly. The fast and slow growth phases of 
IZhao et al. lii^ correspond to growth rates in the matter 
dominated era of M = 6 and M = 1.5, respectively. Our 
result that the core features appear as the growth rate slows 
supports this inside-out view of halo growth. 

We conclude the structural and dynamical features of 
Halo A are valid for haloes near the filtering scale of WDM 
cosmologies in general and form as haloes transition from 
fast growth to the slow accretion phase, remaining stable 
thereafter. 


5 ORIGIN OF THE CORE 

The simulations in this work have shown changes in the 
density and dynamical profiles in the core is common for 
haloes near the filtering mass in WDM cosmologies. The 
previous section determined the formation time of the core 
and demonstrated its stability with time. In this section we 
investigate clues to the core’s origin. 

We examine the cores of Halo A in the extended small 
scale simulations and the cores of Halo B and C at large 
mass scale. We define the proper radius of the core as 60 pc 
for Halo A, 680 pc for Halo B, and 400 pc for Halo C. These 
radii were held constant in proper length and the number of 
particles within each core was calculated for all cosmologies 
at all times a halo progenitor was identified. The top row of 
Figure [ 9 ] shows the growth in core particles in the CDM and 
WDM simulations. It is clear the core forms quickly in WDM 
while the core is built-up more gradually in CDM. The num¬ 
ber of core particles remain approximately constant after 
formation in both CDM and WDM. The epoch of WDM core 
formation occurs shortly after the halo virial mass catches 
up to the CDM halo and the growth rate slows to the CDM 
rate. For example, in the 6 keV small mass scale simulation 
of Halo A the core forms at a ~ 0.075, while Figure [T] shows 
the growth rate slows to approximately the CDM rate at 
about the same time. This agrees with the formation epoch 
of the core feature in the Q profiles. 

We characterize the angular momentum of each core 
with the dimensionless spin parameter, Eqn.[5l and compare 
the evolution of each core spin in the bottom row of Figure[^ 
There is a clear trend between the number of core particles 
and core spin. All cores tend to lose angular momentum over 
time but cores with fewer particles consistently have higher 
spin. There is also a trend that as the filtering scale increases 
the number of core particles first increases then decreases. 
This is most evident in Haloes B and C and is again likely 
due to the later collapse epochs of these haloes. 

In Figure[T0]we show the core particles of Halo A in the 
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I -u.o ■ 

log radius [kpc] 
Halo A-Small-4 keV 



Halo B • Large • CDM 



Halo B • Large -1.1 keV 



Halo B - Large - 0.75 keV 



Halo C • Large • CDM 



Halo C • Large -1.1 keV 
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- 0.75 keV 
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Figure 8 . Evolution of the phase space density profiles of Halo A at small scale and Haloes B and C at large scale (left to right). CDM 
cosmologies are plotted in the top row, WDM in the bottom rows with filtering scale increasing top to bottom. The phase space density 
profiles have been multiplied by to show deviations from power law behaviour. The CDM profiles are consistent with power laws with 
minor fluctuations while the WDM profiles show prominent features in their cores that form shortly after the growth rate slows and are 
stable thereafter. Profiles are plotted with variable colour simply for visual clarity. 


small mass scale simulations at a = 0.05, well before core 
formation. The imag es are centred on t he particles’ centre 
of mass. As noted bv iBusha et al.l ll2007l) . the filtered power 
spectra cause what were multiple clumps in CDM to become 
one collapsing clump in WDM. The core particles addition¬ 
ally become more symmetrically distributed around the cen¬ 
tre of mass as the cosmology becomes warmer. 


This evidence argues for angular momentum playing 
an important role in determining the structure of the core. 
The importance of angular momentum for the shape of the 
inn er profile has bee n emphasized by a number of stud- 
iIHuss et al.lll999'9B lHiotelisl[2002l: lAscasibar et al.|[2004 


ILu et all 200d ). Purely radial orbi ts give a steep inner pro¬ 

file, p oc (lBertschingei|[l98^ . As the amount of angu¬ 

lar momentum is increased particles remain closer to their 
maximum orbital radii resulting in shallower density profiles. 
Angular momentum is dominated by the tangential compo¬ 
nent of the velocity dispersions which are acquired dynami¬ 
cally in both the CDM and WDM simulations since thermal 
velocities were not added to the WDM particles. Interac¬ 
tions with substructure and the global tidal field produce 
tangential components to the particle velocities. An alter¬ 


nativ e possibility is radial orbit instability llBelokurov et ah] 
I2OO8I) . A detailed study of particle orbits is outside the scope 
of this paper, however, after particles collapse the virializa- 
tion process isotropizes their velocities equally well in both 
CDM and WDM as seen from Figure [5] The higher accre¬ 
tion rates in the WDM fast growth phase may al so play a 
role in generating the core as seen in the models of iLu et al.l 

1I2OO6D . 


6 DISCUSSION 

We found the inner structures of dark matter haloes in cos¬ 
mologies with truncated power spectra deviate from their 
profiles in non-truncated cosmologies with mass displaced 
from the intermediate regions to the centre. We have al¬ 
ready shown how our results are in agreement with pre¬ 
vious work on CDM microhaloes Jlshivama et al.l I2OIC1I : 
lAnderhalden fc Diemandll2013l : Ilshivamall20l4) . in this sec¬ 
tion we discuss how our work compares to other WDM and 
HDM studies. 

Investigations of Milky Way satellites in 1—4 keV WDM 
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Scale factor 

Figure 9. Evolution of the number of particles within the core (top) and the spin parameter of the core (bottom) in the extended small 
mass scale simulations of Halo A and large mass scale simulations of Haloes B and C (left to right). The core radius was held fixed in 
proper length at 60 pc for Halo A, 680 pc for Halo B, and 400 pc for Halo C. 



CDM 




4 keV 


6 keV 


Figure 10. Positions at a = 0.05 of core particles in select cos¬ 
mologies of the small mass scale simulations of Halo A. Images 
are centred on the centre of mass of core particles. 


cosmologies have shown the maximum circular velocity de¬ 
creases and the radiu s where this occur s increases for dwarf 
galaxy-sized haloes iLo vell et all I2OI2I : lAnderhalden et 'aD 
I2OI3I : IPolisenskv fc Ricottill201'3 l. Figure [TT] shows the cir¬ 
cular velocity profiles of the small mass scale simulations 
of Halo A at the normalization time. It is clear the rear¬ 
rangement of mass in the inner regions has not affected the 
maximum circular velocity or its location to be in disagree¬ 
ment with the conclusions of previous works. 


It is well established that in WDM cosmologies haloes 
below the truncation scale form later and have lower concen¬ 
trations than CDM halo e s of similar si z e llAvila-Reese et al.l 
2001 L Bode_et_^ 2001: iKnebe et ^ I2OO2I') . In” contrast, 


Anderhalden fc Diemandl ( 2013ll examined their three CDM 


microhaloes and found higher concentrations when NFW 
profiles were fit to the haloes in the filtered power spec¬ 
trum simulations. However, they found fitting generalized 
profiles, where the inner slope is allowed to be a free param¬ 
eter, shifted the scale radius Vs outward and reduced the 
concentrations. 


We investigated the halo concentrations in the low reso¬ 
lution, large refinement volume simulations described in Sec¬ 
tion 2) These simulations were run for the small and large 
box sizes and 6 keV and 1.1 keV cosmologies, respectively. 
We fit NFW and generalized density profiles to the nine 
largest haloes in the refinement volume that were composed 
of > 10® particles and satisfied the relaxation metrics at the 
normalization times. We found seven of the nine haloes had 
greater NFW concentrations in WDM than CDM. However, 
six of the nine halos had lower concentrations when gener¬ 
alized profiles were fit, consistent across box sizes and con¬ 
sistent with the micr ohalo results llAnderhalden fc Pieman^ 
I2OI3I : Iishivamair2014l f. We found the same results when com¬ 
paring the concentrations of Haloes A, B, and C at the end of 
the high-resolution extended simulation set. Interestingly, in 
the simulations with M/ ~ 1.7M the NFW concentrations 
of Haloes B and C were less than their CDM concentrations 
while Halo A was greater. This is consistent with the den¬ 
sity profiles in Figure [3 and the later formation epochs of 
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Haloes B and C, that enhance the filtering of substructure 
for a given filtering scale. 

Our WDM halo prof il es ar e sim ilar to the profile s 
seen by I Avila- Reese et ^ (1200 il l and IColfn et al.l ll2008ll . 
IColfn et al.l I 2008l l simulated five galactic-sized haloes in 
WDM and found profiles were steeper and denser in the 
inner region than the best-fitting NFW profiles. Addition¬ 
ally, our simulati ons explain a feature se en in the HDM clus¬ 
ter simulations of lWang fc Whit3 (l2009l l. The stacked phase 
space density profiles of their 20 most massive haloes show 
a flattening in the inner 0.05Rvir in HDM similar to the 
features seen in the warmest simulations of Figure [HI 

Our re s ults seem to conflict with the work of 
iBusha et al.l (l2007ll . They evolved their CDM and WDM 
simulations into the future until the scale factor a = 100. 
Past the current epoch the cosmological constant quickly 
dominates the density of the universe (Ha —>■ 1) leading 
to exponential expansion and the supression of structure 
growth at a ~ 3. Thus, examining halo properties in the 
far future guarantees the haloes have ample time to relax 
into their equilibrium states. They find the average density 
profiles for haloes near the filtering mass are well fit by the 
NFW profile for r > 0.05i?«ir with only lower concentrations 
below the filtering mass (although the outer slope in both 
CDM an d WDM is ste eper due to the inflating universe, as 
noted bv lRicottilbOOSl l. However, there are two things that 
complicate comparison of these simulations to ours. First is 
the difficulty due to the different epochs the haloes are ex¬ 
amined at. We examine our haloes shortly after the end of 
the fast accretion phase when th e inner profile is set but the 
outskirts are still growing while iBusha et al.l examine their 
haloes well after all halo growth has stopped and R^ir has 
reached a maximum. The effects we see in the inner halo will 
therefore be at radii sma ller than the convergence radius in 
the haloes simulated by iBusha et al.l . Also, they explored 
haloes at a larger mass scale and used a much greater filter¬ 
ing mass, 1.2 x IO^^Mq. In this work we have shown how 
even the density profiles of CDM haloes have a dependence 
on mass with smaller differences between CDM and WDM 
profiles for larger halo masses, for a fixed ratio of the filtering 
mass to halo mass. 

Finally, we comment on the effects of adding ther¬ 
mal velocities appropriate to the adopted WDM models to 
the simulation particles. Thermal WDM particles decou¬ 
ple with a finite fine-grain phase space density that im¬ 
poses an upper limit on their density, resulting in soft cores 
in collapsed haloes. The radius of this core depends on 
the mass of the WDM pa rticle and the mass of the halo 
(iHogan fc Dalcantonll200(]l l. For the warmest cosmologies of 
Halo A the core radius is ~ 4 x 10~^R vir which agrees wit h 
the core sizes seen in the simulations of iMaccio et al.l ll2012ll . 
The thermal core would be about the size of the adopted 
softening lengths, far below the scales where the WDM pro¬ 
files deviate from the CDM profile. 


7 SUMMARY 

We tested the claim that the virialization process erases 
all information about the initial conditions and produces 
universal profiles in gravitationally collapsed dark matter 
haloes. We simulated an isolated halo with an early for- 



Figure 11. Circular velocity profiles of the small mass scale sim¬ 
ulations of Halo A. The WDM profiles are grouped and plotted 
against the CDM profile for clarity. 


mation epoch at three mass scales from 2 x IO^Mq to 
6 X 1O^°M0 in CDM and a variety of WDM cosmologies 
where the formation of structures below the hltering scale 
is suppressed. We examined when the halo was composed of 
~ 10^ particles at each scale and the halo was in the slow 
growth phase. We studied two additional haloes to account 
for cosmic variance. We found the haloes were changed both 
structurally and dynamically by the truncation in the WDM 
power spectra. Substructures in the mass range 
have a detectable effect on the slope of the inner density pro¬ 
file in haloes because when substructure on this 

mass scale is filtered out in the 34 keV WDM cosmology the 
profile visibly deviates from CDM, becoming steeper. We 
summarize our main findings below. 

• Density profiles are steeper in WDM than CDM for 
haloes near the filtering scale. For Mf < M, the slopes ap¬ 
proach a constant value moving from the virial radius to the 
inner h alo, whose value is well described by the htting for¬ 
mula of llshivamal (|2ni4l l. The value of the constant slope and 
its radial span depend on the filtering scale, moving outward 
and becoming steeper as the ratio of the filtering to halo 
mass gets larger. Our simulations also examined the regime 
Mf > M and found the density profiles do not approach a 
constant asymptotic slope at intermediate radii but continue 
to flatten monotonically. At smaller radii, for all the masses 
and WDM cosmologies, the slopes do not remain at the con¬ 
stant value found at intermediate radii but grow shallower 
and become softer than CDM for r < 0.01R„ir, in agreeme nt 
with halo models llWilliams et al.|[2004l : IVifias et al.ll2012l ~). 

• Particle velocity dispersions increase in the inner pro¬ 
files while velocity anisotropies after virialization are largely 
similar across cosmologies. The changes in density and veloc¬ 
ity dispersion create a deviation from a simple power law in 
the inner phase space density profiles. This deviation reaches 
a maximum at a radius proportional to the filtering scale. 

• The core features of the profiles are set early when the 
halo mass growth rate M ~ 6 — 10. At higher growth rates 
the profiles exhibit large fluctuations. After formation the 
core features are dynamically stable. 

• The halo mass structure is rearranged in WDM com¬ 
pared to CDM with radii < 0.1R„ir gaining mass at the 










































expense of radii 0.1 — 0.47?„ir. Furthermore, the core is built 
up gradually in CDM from particles distributed asymmet¬ 
rically in clumps about the centre of mass in contrast to 
WDM where the core is formed in an impulsive event from 
particles distributed smoothly and symmetrically. There is 
a correlation between the number of particles and the spin 
of the core with more massive cores having less spin. This 
argues for angular momentum as the physical mechanism re- 
sp onsible for the difference s in profiles as seen in the models 
of iDel Popol3 (l2009l. 120101) . 

• We found a dependence on mass in the CDM profiles 
with la r ger ha loes exhibiting a steeper density profile as in 
iRicottil (I2OO3I ). The spin parameter decr eases with increas - 
ing mass in agreement with the models of lPel Popolol (l2009l ) 
that more massive haloes have less angular momentum re¬ 
sulting in steeper profiles. However, the effects are much 
weaker than those of the truncated power spectrum. The 
WDM haloes had similar spins across mass scales and also 
had similar profiles. 


Our work complements and reinforces the results of 
recent simulations of microhaloes near the CDM filtering 
scale (lls hivama et al.l[201ol : lAnderhalden fc Diemandll2013l : 
llshivama 2014 ) and shows that the shape of halo profiles 
cannot be parameterized simply by a generalized NFW or 
Einasto profile with a concentration or scale radius de¬ 
pendent on the mass or cosmology. The halo shape is 
more complex, with logarithmic slopes that can vary non- 
monotonically and with features in the profile that reflect 
the shape of the initial matter power spectrum. Thus, in 
general haloes cannot be fitted by a universal density profile. 
This is actually good news because it may become feasibile 
to find fingerprints of the initial power spectrum of perturba¬ 
tions on galactic or sub-galactic scales in the density profiles 
of dark matter dominated dwarf galaxies or clusters. 


ACKNOWLEDGMENTS 

The simulations presented in this work were run on the 
DEEPTHOUGHT computing cluster at the University of 
Maryland College Park, the Cray XE6 GARNET at the U.S. 
Army Engineer Research and Development Center and the 
SGI Ice X SPIRIT at the U.S. Air Force Research Labo¬ 
ratory. Basic research in astrophysics at NRL is funded by 
the U.S. Office of Naval Research. EP acknowledges support 
under the Edison Memorial Graduate Training Program 
at the Naval Research Laboratory. The authors acknowl¬ 
edge the University of Maryland supercomputing resources 
(http://www.it.umd.edu/hpcc) made available in conduct¬ 
ing the research reported in this paper. MR’s thanks the 
National Science Foundation and NASA for support under 
the grants NASA NNXIOAHIOG, NSF CMMI1125285 and 
the Theoretical and Computational Astrophysics Network 
(TCAN) grant AST1333514. 


REFERENCES 

Anderhalden D., Diemand J., 2013, JCAP, 4, 9 
Anderhalden D., Schneider A., Maccib A. V., Diemand J., 
Bertone G., 2013, JGAP, 3, 14 


Fingerprints of the ICs 13 

Ascasibar Y., Yepes G., Gottlober S., Muller V., 2004, MN- 
RAS, 352, 1109 

Austin G. G., Williams L. L. R., Barnes E. I., Babul A., 
Dalcanton J. J., 2005, ApJ, 634, 756 
Avila-Reese V., Colfn P., Valenzuela O., D’Onghia E., Fir- 
mani G., 2001, ApJ, 559, 516 
Barnes E. L, Williams L. L. R., Babul A., Dalcanton J. J., 
2006, ApJ, 643, 797 

Belokurov V., Walker M. G., Evans N. W., Faria D. G., 
Gilmore G., Irwin M. J., Koposov S., Mateo M., Olszewski 
E., Zucker D. B., 2008, ApJ, 686, L83 
Bertschinger E., 1985, ApJS, 58, 39 
Bertschinger E., 2001, ApJS, 137, 1 
Bode P., Ostriker J. P., Turok N., 2001, ApJ, 556, 93 
Bullock J. S., Kravtsov A. V., Weinberg D. H., 2001, ApJ, 
548, 33 

Busha M. T., Evrard A. E., Adams F. C., 2007, ApJ, 665, 
1 

Golfn P., Avila-Reese V., Valenzuela O., 2000, ApJ, 542, 
622 

Golfn P., Valenzuela O., Avila-Reese V., 2008, ApJ, 673, 
203 

Dekel A., Devor J., Hetzroni G., 2003, MNRAS, 341, 326 
Del Popolo A., 2009, ApJ, 698, 2093 
Del Popolo A., 2010, MNRAS, 408, 1808 
Diemand J., Moore B., Stadel J., 2004, MNRAS, 353, 624 
Einasto J., 1965, Trudy Astrofizicheskogo Instituta Alma- 
Ata, 5, 87 

Eisenstein D. J., Hu W., 1998, ApJ, 496, 605 
Gott HI J. R., 1975, ApJ, 201, 296 

Graham A. W., Merritt D., Moore B., Diemand J., Terzic 
B., 2006, AJ, 132, 2701 
Gunn J. E., Gott III J. R., 1972, ApJ, 176, 1 
Hiotelis N., 2002, A&A, 382, 84 

Hogan G. J., Dalcanton J. J., 2000, Phys. Rev. D, 62, 
063511 

Huss A., Jain B., Steinmetz M., 1999a, ApJ, 517, 64 
Huss A., Jain B., Steinmetz M., 1999b, MNRAS, 308, 1011 
Ishiyama T., 2014, ApJ, 788, 27 

IshiyamaT., Makino J., Ebisuzaki T., 2010, ApJ, 723, L195 
Jing Y. P., 2000, ApJ, 535, 30 
Jing Y. P., Suto Y., 2000, ApJ, 529, L69 
Klypin A. A., Trujillo-Gomez S., Primack J., 2011, ApJ, 
740, 102 

Knebe A., Devriendt J. E. G., Mahmood A., Silk J., 2002, 
MNRAS, 329, 813 

Knollmann S. R., Knebe A., 2009, ApJS, 182, 608 
Lacey C., Cole S., 1993, MNRAS, 262, 627 
Lovell M. R., Eke V., Frenk C. S., Gao L., Jenkins A., The- 
uns T., Wang J., White S. D. M., Boyarsky A., Ruchayskiy 

O. , 2012, MNRAS, 420, 2318 

Lu Y., Mo H. J., Katz N., Weinberg M. D., 2006, MNRAS, 
368, 1931 

Ludlow A. D., Navarro J. F., Angulo R. E., Boylan-Kolchin 
M., Springel V., Frenk G., White S. D. M., 2014, MNRAS, 
441, 378 

Ludlow A. D., Navarro J. F., Boylan-Kolchin M., Bett 

P. E., Angulo R. E., Li M., White S. D. M., Frenk G., 
Springel V., 2013, MNRAS, 432, 1103 

Ludlow A. D., Navarro J. F., Li M., Angulo R. E., Boylan- 
Kolchin M., Bett P. E., 2012, MNRAS, 427, 1322 
Lynden-Bell D., 1967, MNRAS, 136, 101 
























14 E. Polisensky and M. Ricotti 


Maccio A. V., Dutton A. A., van den Bosch F. C., 2008, 
MNRAS, 391, 1940 

Maccio A. V., Dutton A. A., van den Bosch F. C., Moore 

B. , Potter D., Stadel J., 2007, MNRAS, 378, 55 
Maccio A. V., Paduroiu S., Anderhalden D., Schneider A., 

Moore B., 2012, MNRAS, 424, 1105 
Manrique A., Raig A., Salvador-Sole E., Sanchis T., Solanes 
J. M., 2003, ApJ, 593, 26 

Moore B., Ghigna S., Governato F., Lake G., Quinn T., 
Stadel J., Tozzi P., 1999, ApJ, 524, L19 
Moore B., Quinn T., Governato F., Stadel J., Lake G., 
1999, MNRAS, 310, 1147 

Navarro J. F., Frenk C. S., White S. D. M., 1996, ApJ, 462, 
563 

Navarro J. F., Frenk C. S., White S. D. M., 1997, ApJ, 490, 
493 

Navarro J. F., Hayashi E., Power G., Jenkins A. R., Frenk 
G. S., White S. D. M., Springel V., Stadel J., Quinn T. R., 
2004, MNRAS, 349, 1039 

Neto A. F., Gao L., Belt P., Cole S., Navarro J. F., Frenk 

C. S., White S. D. M., Springel V., Jenkins A., 2007, MN¬ 
RAS, 381, 1450 

Nusser A., Sheth R. K., 1999, MNRAS, 303, 685 
Peebles P. J. E., 1974, ApJ, 189, L51 
Polisensky E., Ricotti M., 2011, Phys. Rev. D, 83, 043506 
Polisensky E., Ricotti M., 2014, MNRAS, 437, 2922 
Power C., Navarro J. F., Jenkins A., Frenk C. S., White 
S. D. M., Springel V., Stadel J., Quinn T., 2003, MNRAS, 

338, 14 

Prada F., Klypin A. A., Cuesta A. J., Betancort-Rijo J. E., 
Primack J., 2012, MNRAS, 423, 3018 
Press W. H., Schechter P., 1974, ApJ, 187, 425 
Ricotti M., 2003, MNRAS, 344, 1237 
Ricotti M., Pontzen A., Viel M., 2007, ApJ, 663, L53 
Ricotti M., Wilkinson M. L, 2004, MNRAS, 353, 867 
Sommer-Larsen J., Dolgov A., 2001, ApJ, 551, 608 
Springel V., 2005, MNRAS, 364, 1105 
Subramanian K., Gen R., Ostriker J. P., 2000, ApJ, 538, 
528 

Syer D., White S. D. M., 1998, MNRAS, 293, 337 
Taylor J. E., Navarro J. F., 2001, ApJ, 563, 483 
Vinas J., Salvador-Sole E., Manrique A., 2012, MNRAS, 
424, L6 

Wang J., White S. D. M., 2007, MNRAS, 380, 93 
Wang J., White S. D. M., 2009, MNRAS, 396, 709 
Wechsler R. H., Bullock J. S., Primack J. R., Kravtsov 
A. V., Dekel A., 2002, ApJ, 568, 52 
White S. D. M., 1996, in Lahav O., Terlevich E., Terlevich 
R. J., eds. Gravitational dynamics Violent Relaxation in 
Hierarchical Clustering, p. 121 
Williams L. L. R., Babul A., Dalcanton J. J., 2004, ApJ, 
604, 18 

Zhao D. H., Jing Y. P., Mo H. J., Borner G., 2003, ApJ, 
597, L9 

Zhao D. H., Jing Y. P., Mo H. J., Borner G., 2009, ApJ, 
707, 354 

Zhao D. H., Mo H. J., Jing Y. P., Borner G., 2003, MNRAS, 

339, 12 


